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We present preliminary results from the QCDSF/UKQCD collaborations for the second moment of the pion's 
distribution amplitude with two flavours of dynamical fermions. We use nonperturbatively determined renormal- 
isation coefficients to convert our results to the MS scheme at 5 GeV 2 . Employing a linear chiral extrapolation 
from our large pion masses > 550 MeV, we find (£ 2 ) = 0.281(28), leading to a value of 02 = 0.236(82) for the 
second Gegenbauer moment. 



1. INTRODUCTION 

The distribution amplitude of the pion, 
4>n(x, fj?), contains information on how the pion's 
longitudinal momentum is divided between its 
quark and anti-quark constituents when probed 
in diffractive dijet production at E-791 [I] and 
exclusive pion photoproduction at CLEO 2 . 

The pion's distribution amplitude is defined on 
the light cone as 

(0\d(0) 7fll5 u(z)\ir + (p)) = 

ifnPn / dx e- lxp ' z (t>^(x,n) , (1) 
Jo 

where z 2 = is a vector along the light cone, x is 
the fraction of the pion's longitudinal momentum, 
p, carried by the quark (x = 1 — x for the anti- 
quark), / w is the pion decay constant, and /1 is 
the factorisation scale. On a Euclidean lattice, 
we are not able to compute matrix elements of 
bilocal operators such as d(0)j tl j 5 u(z) in Eq. JIJ, 
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instead we make use of the light-cone operator 
product expansion which allows one to calculate 
Mellin moments of Eq. JQ) via the computation 
of matrix elements of local operators. The n th 
moment of the pion's distribution amplitude is 
defined as 

(C) = J d£f l <M£,Q 2 ) , (2) 

where £ = x — x = 2x — 1, and can be extracted 
from the matrix elements of twist-2 operators 

<0|O {M0 ...„ B} (0)|7r(p)> = U Pll0 . . .Pp„(C> + - • • , 

(3) 

where 

O lto ... ltn (0) = H)"?7«,75 B MI • B M „ V- , (4) 

D=D — D and {• • • } denotes symmetrisation of 
indices and subtraction of traces. We implement 
the standard normalisation by setting = 1. 
Due to G-parity the first moment, vanishes 
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for the pion, hence the first nontrivial moment 
that we are able to calculate is (£ 2 ). 

Although the first lattice calculations of (£ 2 ) 
appeared almost 20 years ago there has been 
surprisingly little activity in this area in recent 
times j'll5l6| to complement other theoretical in- 
vestigations, e.g. |7l8l9l1 1)11111 '211411 411 511 fill VI 
f 8 . The current state-of-the-art lattice calcu- 
lation comes from Del Debbio et al. |S| who 
performed a simulation in quenched QCD and 
renormalised their results perturbatively to the 
MS scheme at /j = 2.67 GeV, (£ 2 ) W = 
2.67 GeV) = 0.280(49)^;^ . 

In these proceedings we present preliminary re- 
sults from the QCDSF/UKQCD collaborations 
for (£ 2 } MS in two flavour lattice QCD. These re- 
sults complement our preliminary results on the 
pion form factor |19| . 
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Table 1 

Lattice parameters: Gauge coupling (3, sea quark 
hopping parameter K SC a, lattice volume, lattice 
spacing and pion mass. 



2. OPERATORS 

The -ff (4) representation on the lattice leads us 
to consider two operators which we call O a and 
O b [HI, e.g. 



Eq. JSJ can mix with the following operator [3U| 



°{4i2} = -J d {t d i (^72}75^) 



(7) 



P = (p,p,o)- 

^412 - {412} , (5) 

P = (p,0,0): 

nb (n ^{422} + 0(433} \ 
0411 = I 0{411} ^ ) ' ( b > 

where p — 2ir/L s and L s is the spatial extent of 
our lattice. From Eq. J3J), we see that O a requires 
two spatial components of momentum, while O b 
needs only one. Consideration of this fact alone 
would lead one to choose O b , since units of mo- 
menta in different directions on the lattice lead to 
a poorer signal. However, lattice operators with 
two or more covariant derivatives can mix with 
operators of the same or lower dimension. For 
forward matrix elements, O b suffers from such 
mixings while O a does not. For matrix elements 
involving a momentum transfer between the two 
states, i.e. nonforward matrix elements, both op- 
erators O a and O b can mix with operators involv- 
ing external ordinary derivatives, i.e. operators of 



the form d^dy ■ 



(tp ■ ■ ■ ip). For example, 0% 
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The situation for O b is a lot worse as it can po- 
tentially mix with up to seven different operators 
|2*D|. Hence a complete calculation of (£ 2 ) would 
require knowledge of the mixing coefficients and 
renormalisation constants for all of these mixing 
operators. It now becomes obvious that O a offers 
the best possibility to extract a value of (£ 2 ) from 
a lattice simulation. 

Although the mixing coefficient for O a ' dd is not 
yet known, we expect that it is small and hence 
we anticipate that the contribution to (£ 2 ) from 
Qa, dd w jvj k e srna rj Hence, for the rest of the 
work presented here, we will consider only the 
contribution from O a . In a forthcoming publica- 
tion, we will attempt to address all mixing issues 
associated with both operators O a and O b . 

3. LATTICE TECHNIQUES 

We simulate with Nf = 2 dynamical configu- 
rations generated with Wilson glue and nonper- 
turbatively 0{a) improved Wilson fermions. For 
four different values f3 = 5.20, 5.25, 5.29, 5.40 
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Figure 1. Ratio defined in Eq. using O a from 
Eq. J3|j, on a 24 3 x 48 lattice with = 5.40 and 

^sea — 

0.13560. 
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Figure 2. (£ 2 ) as a function of m\ for O a from 
Eq. 0). 



and up to four different kappa values per beta we 
have generated 0(2000 — 8000) trajectories. Lat- 
tice spacings and spatial volumes vary between 
0.075-0.123 fm and (1.5-2.2 fm) 3 respectively. A 
summary of the parameter space spanned by our 
dynamical configurations can be found in Tabled 
We set the scale via the force parameter, with 
ro = 0.5 fm. 

Correlation functions are calculated on config- 
urations taken at a distance of 10 trajectories us- 
ing 4 different locations of the fermion source. We 
use binning to obtain an effective distance of 20 
trajectories. The size of the bins has little ef- 
fect on the error, which indicates residual auto- 
correlations are small. 

We calculate the average of matrix elements 
computed with three choices of pion momenta 
Po = (p,P,0), Pi = (p,0,p), P2 = (0,P,p), with 
the indices of the operators (Eq. (JSJ) chosen ac- 
cordingly. 

We define a pion two-point correlation function 

as 

C°(t,p) = ^ e ^(o(a?,t)J(0,0)t) , 

X 

- ^(0\O(0Mp))e- Et , t»0, (8) 

where Z — (n(p)\ J(0)t|0) and we use J{x) = 
7r(x) = il)(x)^5%l)(x) as the interpolating opera- 



tor for the pion. The matrix elements in Eq. 
are then extracted from the following ratios of 
two-point functions, 

where i and j are spatial indices, and O4 is the 
operator given in Eq. @ with no derivatives and 
Mo = 4. 

Figure ^ shows a typical example of the ratio 
in Eq. © where we clearly observe two plateaus 
for t < Lt/2 and t > Lt/2, where Lt is the time 
extent of the lattice. After extracting R a from 
the plateaus, we use Eq. © to extract (£ 2 ). 

In general, bare lattice operators must be 
renormalised in some scheme, S, and at a scale, 
M, 

O s {M) = Z$(M)O bare , (11) 

so in order to calculate a renormalised value for 
(£ 2 ); we must consider 

(ef(M) = ^MMl { e) hare . (12) 

We choose to renormalise to the MS scheme at 
a scale of /i 2 = 5 GeV 2 . Further details of our 
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Figure 3. Chiral extrapolation at constant (3 for 
(3 = 5.40 (top) and (3 = 5.29 (bottom) for O a 
from Eq. 0). 



Figure 4. Chiral extrapolation at constant (3 for 
(3 = 5.25 (top) and (3 = 5.20 (bottom) for O a 
from Eq. 0). 



renormalisation techniques can be found in |21| 
and a forthcoming publication. 

4. RESULTS FOR (£ 2 ) 

For each of our datasets, we extract a value 
for (£ 2 }bare from Eq. © and renormalise using 
Eq. I|12|l . Figure |21 shows these results plotted as 
a function of m 2 . Here we observe that the results 
are approximately constant as we vary the pion 
mass. 

In order to obtain a result in the continuum 
and chiral limits, we first extrapolate our results 
at constant (3 to the chiral limit. In Fig.|2]we dis- 
play the chiral extrapolations for (3 = 5.40 (top) 
and 5.29 (bottom), while Fig. 0] contains the cor- 
responding extrapolations for (3 — 5.25 (top) and 
5.20 (bottom). These results exhibit only a mild 
dependence on the quark mass and their values 
in the chiral limit agree within errors. 

Now that we have calculated results in the chi- 
ral limit for each choice of /3, we are in a position 
to examine the behaviour of our results as a func- 
tion of the lattice spacing. In Fig. El we use the 
values of ro calculated in the chiral limit for each 
13 (see Table 3 of Ref. [22]) to study the depen- 
dence of our results on the lattice spacing. Here 



we observe that even though our operators are 
not O(a)-improved, we find a negligible depen- 
dence on the lattice spacing, at least when com- 
pared to the statistical errors. Hence we take the 
result at our smallest lattice spacing (largest (3) 
as our result in the continuum limit. 

We find in the continuum limit at the physi- 
cal pion mass, the second moment of the pion's 
distribution amplitude to be 

(?) m (v 2 = 5 GeV 2 ) = 0.281(28) , (13) 

which is very close to the value (£ 2 ) M (/i = 
2.67 GeV) = 0.280(49)1^?° found in Ref. 0, 
and larger than the asymptotic value, /i 2 — > oo 

0»(O = |(l-e 3 ) => <a = 0.2. (14) 

5. GEGENBAUER MOMENT a 2 

The distribution amplitude, </>(:r,/! 2 ), can be 
expanded in a series of even Gegenbauer poly- 
nomials, C| n (2x - 1) [EE], 

OO 

0(x, M 2 ) = 6x(l-x)^a 2 „( M 2 )c|„(2 2 :-l), (15) 
n=0 
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Figure 5. Results for each value of (3 in the chiral 
limit as a function of a 2 for O a from Eq. f3|) . 

where a 2n (/it 2 ) are the multiplicatively renormal- 
isable Gegenbauer moments. Since the higher 
moments, n > 4, are expected to be small, it 
is common to truncate Eq. (|15|) 

<j)(x,fi 2 ) = 6x(l - x) jl + a 2 {fi 2 )cl (2x - 1) 

+ a 4 (n 2 )cl (2x - 1)} . (16) 

In the asymptotic limit, [i 2 — > oo, all a 2 „(/i 2 ) = 0, 
for n > 0. 

Analyses of CLEO data for F nil * constrain the 
pion distribution amplitude by calculating a re- 
lationship between the first two Gegenbauer mo- 
ments (see eg., I15I16) 1 ). together with upper and 
lower bounds on their respective values. 

Taking the second £ = 2x — 1 moment of the 
r.h.s. of Eq. Ijl6|l gives 

<£ 2 } - f d£ M 2 ) = ^ (7 + 12a 2 ( M 2 )) , 

(17) 

which allows us to extract a 2 , but not 04. In order 
to place a constraint on 04, we would need to cal- 
culate the fourth moment of the pion's distribu- 
tion amplitude which requires using an operator 
involving four covariant derivatives. 

Using our result in Eq. I|13|l we calculate 

a 2 {fi 2 = 5 GcV 2 ) = 0.236(82) , (18) 



which is larger than the values a 2 (lGeV) = 
0.07(1) [HI and a 2 (lGeV) = 0.19(19) [TBI, but 
within the usual constraints < a 2 (l GeV) < 0.3 
commonly quoted in the literature, even when we 
consider that our value will increase slightly when 
we run our result to a smaller scale of y, = 1 GeV. 

6. CONCLUSIONS 

We have presented a preliminary result for the 
second moment of the pion's distribution am- 
plitude calculated on lattices generated by the 
QCDSF/UKQCD collaboration with two flavours 
of dynamical fermions. We use nonperturbatively 
determined renormalisation coefficients to con- 
vert our result to the MS scheme at 5 GeV 2 . We 
find (£ 2 ) = 0.281(28), which is very close to an 
earlier lattice result and larger than the asymp- 
totic value. 

Using a fourth order Gegenbauer polynomial 
expansion, we calculate a value for the sec- 
ond Gegenbauer moment, a 2 (^ 2 = 5 GeV 2 ) = 
0.236(82). 

Although we have only employed a linear chi- 
ral extrapolation and our operators are not 0(a)- 
improvcd, the chiral and continuum extrapola- 
tions do not seem to be a major source of sys- 
tematic error when compared to the statistical 
errors. These issues will be addressed in more de- 
tail in a forthcoming coming publication, where 
we also intend to investigate finite size and (par- 
tially) quenching effects as well as renormalisation 
group running of the relevant matrix elements. 
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